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ABSTRACT 

We explore the growth of super-massive black holes and host galaxy bulges in the 
galaxy population using the Millennium Run ACDM simulation coupled with a model 
of galaxy formation. We find that, if galaxy mergers are the primary drivers for both 
bulge and black hole growth, then in the simplest picture one should expect the ms h — 
mbuige relation to evolve with redshift, with a larger black hole mass associated with a 
given bulge mass at earlier times relative to the present day. This result is independent 
of an evolving cold gas fraction in the galaxy population. The evolution arises from 
the disruption of galactic disks during mergers that make a larger fractional mass 
contribution to bulges at low redshift than at earlier epochs. There is no comparable 
growth mode for the black hole population. Thus, this effect produces evolution in the 
mBH - rribuige relation that is driven by bulge mass growth and not by black holes. 
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1 INTRODUCTION 

Super-massive black hole masses are strongly correlated 
with their host bulge stellar mass, th e so-called ttibh - 
mbuige relation iMagorrian et ahl Il99fj : iMarconi fc Hu nil 
l2003r lHaring fa T^fHoOdl) . This*ls at least true in the lo¬ 
cal universe, but also expected to extend out to higher red- 
shifts. This correlation suggests a common mechanism link¬ 
ing the growth of these two galactic components, with ev¬ 
idence proposing galaxy mergers as the most likely candi¬ 
date. If true, and given that the global galaxy merger rate 
in a ACDM universe evolves strongly with time, one may 
ask if we should expect to see the tobh - mb u i ge relation also 
evolve. 

Support for the idea that bulges and black holes grow 
through mergers arises primarily from the success of nu¬ 
merical simulations and galaxy formation models in repro¬ 
ducing many observed galaxy scaling relations. Such works 
illustrate that much of the bulge mass of a galaxy can 
be accounted for by the disruption of disk stars from the 
merger proge nitors, and merger triggered starbursts in the 
cold gas d isk ll Barne 3 11992 ); IMihos fc Hernauistlll994l . tl996 : 
ICox et all 12004) 7 As a growth mechanism for black holes, 
merger induced perturbations of the gas close to the cen¬ 
tral massive object can drive gas inward, fueling what is 
obs erved to be a ‘quasar’ period in a galaxy’s history ( see 
e.g. iKauffmann fc HaehnellJl[2000l : iDi Matteo et al.1l200Fiil . 

In this simple picture the amount of cold gas present in 
a merging system plays a large part in how rapidly the black 
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hole and bulge can grow. If the growth dependence for each 
is a simple constant scaling with gas mass, as is commonly 
assumed in many models of galaxy formation, then their 
mass ratio will, on average, be approximately independent 
of any evolution in the global cold gas fraction. This is be¬ 
cause both bulges and black holes then co-evolve at a similar 
pace (drawing their new mass from the same gas reservoir). 
Furthermore, during a merger bulges will add to bulges, and 
black holes may coalesce. Thus, from this alone, one should 
expect little evolution in the rriBH - mbuige relation. 

In this paper we explore an additional growth channel 
through which bulges gain mass that black holes do not have. 
This is the disruption of merged satellite disks, and in the 
event of a major merger, the disruption of the central galaxy 
disk. The stellar mass in such disks will have previously never 
contributed to the uibh -m-buige relation. If the bulge growth 
rate from such disrupted disks is not constant with time, 
then evolution in the tobh - nibuige relation can occur. 

We investigate this behavior using th e Millennium Run 
ACDM simulationjjSgnnge^tahJj^OO^I) and a model for 
galaxy formation lICrotone^Ln^OOfill . This model follows 
the growth of galaxies (including their individual disk, bulge 
and black hole components) from high redshift to the present 
day, and provides a solid framework within which to un¬ 
dertake our analysis. The results we find, however, are not 
be unique to our particular implementation of the galaxy 
formation physics but arise from the simple assumptions 
described above regarding black hole and bulge growth in 
galaxies. Our aim in using this particular model is to illus- 
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trate what one may expect to see if these underlying growth 
mechanisms turn out to be true. 

This paper is organized as follows. In Section |2] we 
briefly describe the Millennium Run ACDM dark matter 
simulation and our model of galaxy formation, including the 
simple implementation of bulge and black hole growth. In 
Section [3 we will use this model to investigate how black 
hole and bulges co-evolve together in the galaxy population 
from high redshift to the present. We finish in Section0]with 
a discussion of the ?tibh -uibuige relation in light of these re¬ 
sults. 


2 GALAXY FORMATION IN A 

COSMOLOGICAL CONTEXT 

The galaxy formation model we use to study th e mbh - 
mbu ige relation is identical to that described in lCroton et alJ 
(‘ 2006 1 (including parameter choices), with the exception 
of one non-essential detail, discussed below. This model of 
galaxy formation is implemented on t op o f the Millennium 
Run ACDM dark matter simulation llSoringel et, al.l feoQH'). 
Below we briefly outline the relevant aspects of the simula¬ 
tion and model to our current work, and refer the interested 
reader to the above references for further information. 

The Millennium Run follows the dynamical evolution of 
10 10 dark matter particles in a periodic box of side-length 
500h -1 Mpc with a mass resolution per particle of 8.6 x 
10 8 /i^'Mg. We adopt cosmological parameter values consis¬ 
tent w ith a combined analysis of the 2dFGRSj|Collesset^ffi| 
|20Wj)_and first year WMAP data llSpergel^taiTl2nO? 
ISehak et al ] 120051) : n A = 0.75, On = f2 dm + n b = 0.25, 
fib = 0.045, h = 0.73, and as = 0.9. Friends-of-friends 
(FOF) halos are identified in the simulation using a linking 
length of 0.2 the mean particle separation, while substruc¬ 
ture within each FOF halo is found with an im proved a nd ex¬ 
tended version of the SUBFIND algorithm of lSpringel et al] 
i200ll ). Having determined all halos and subhalos at all out¬ 
put snapshots we then build the hierarchical merging trees 
that describe in detail how structures grow as the universe 
evolves. These trees form the backbone onto which we couple 
our model of galaxy formation. 

Inside each tree, virialised dark matter halos at each 
redshift are assumed to attract ambient gas from the sur¬ 
rounding medium, from which galaxies form and evolve. Our 
model effectively tracks a wide range of galaxy formation 
physics in each halo, including reionization of the inter- 
galactic medium at high redshift, radiative cooling of hot 
gas and the formation of cooling flows, star formation in the 
cold disk and the resulting supernova feedback, black hole 
growth and active galactic nuclei (AGN) feedback through 
the ‘quasar’ and ‘radio’ epochs of AGN evolution, metal 
enrichment of the inter-galactic and intra-cluster medium, 
and galaxy morphology shaped through mergers and merger 
induced starbursts. As galaxy mergers and the resulting 
growth of bulges and black holes are central to the ques¬ 
tions at hand, we will now describe these in more detail. 


2.1 A simple picture of black hole and bulge 
growth 

2.1.1 The “static” model 

A satellite galax y orbiting within a larger h alo will feel dy¬ 
namical friction llBinnev fc Tremainelll98 7) and eventually 
spiral inward to merge with the central galaxy of the system. 
Mergers are believed to trigger galactic starbursts, where 
some (perhaps large) fraction of the cold disk gas is con¬ 
verted into stars on a timescale much shorter than that typ¬ 
ically found in quiescent star forming disks. To model this 
event, when a merger occurs we assume the following mass of 
stars are formed in a burst from the combined cold gas mass 
of the progenitor g alaxies, as found in the SPH simulations 
of lGox et alJ ( 20041) : 

Am s tarburst — 0.56 777. R 777 C old , (1) 


where w.r = m sa .t/rn ce ntrai is the merger mass ratio of 
the merging galaxies, and m co id the total mass of cold gas 
present during the merger. These stars contribute to the 
spheroid of the final galaxy. For the results presented in this 
paper, the typical mass of stars in a bulge formed through 
starbursts is ~ 10%, which, from Eq. Q indicates an aver¬ 
age gas fraction in the merging progenitors of ^,30%. This 
is consistent with smoothed particle hydrodynamic simu¬ 
lations of merging galaxies which suggest that such a gas 
fract ion is required to explain the local Fundame ntal Plane 
(e.g. iHernauist et al]|l99 3i: iRobei'tsoiT et al . 2005). 

Mergers also perturb the cold gas disk, and this can 
trigger the accretion of g as onto the central super-massive 
black hole. ICroton et aLl f 2006 ) showed that, under reason¬ 
able assumptions, merger triggered ‘quasar’ mode events are 
sufficient to reproduce the local rriBH - mb„i ge relation as well 
as the observed local black hole mass density of the universe. 
To include such events, we app l y an e mpiri cal r elation sim¬ 
ilar to that described in iKauffmann^z Haehnelt ( 12000 ) and 
assume that during a merger the gas accreted onto the black 
hole is proportional to the cold gas present, but in a way that 
is less efficient in lower mass halos: 


Am.BH = 0.03 mn 


1 + (280 kins 1 /V v ir) 2 


Ulcold ? 


( 2 ) 


where VW is the virial velocity of the system, and ttir and 
m co id are defined above. Here, the coefficient 0.03 normal¬ 
izes the m.BH ~ mbuige relation to match that observed lo¬ 
cally. It is important to note that in this picture the ra¬ 
tio AmBH/Am a tarburst, he. the relative growth rate of black 
holes and bulges due to converted cold gas, is expected to be 
essentially constant and independent of redshift, even if the 
gas fraction itself changes with redshift (note that the virial 
velocity of a system is only weakly dependent on time). 

In addition to starbursts, bulges also grow from the 
stellar remnants of merged satellites. In our implementation 
any existing satellite disk is permanently disrupted during 
a merger and its stars, whose orbits will be heavily random¬ 
ized from strong tidal forces, are added to the final galaxy 
bulge. Furthermore, we assume that if the baryonic mass 
ratio of the merging galaxies is large enough a major merger 
has occurred. Major mergers are sufficiently energetic that 
the disk of the central galaxy is also destroyed and its stars 
added to the bulge: we trigger such events when mR>0.3. 
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To summarize, aside from the benign contribution dur¬ 
ing mergers of bulges to bulges and black holes to black 
holes, bulges in our model grow through both starbursts 
and disrupted disks, whereas black holes grow only by ac¬ 
cretion. Importantly, black holes have no comparable growth 
mode from disrupted disks, and in Section 0 we will explic¬ 
itly show the significance of this effect. Finally, given that 
we have made no implicit assumption regarding evolution in 
the growth of either bulges or black holes, we hereafter refer 
to this model as the static model. 


2.1.2 The “dynamic” model 

The above model is not unique in its ability to reproduce the 
local black hole and bulge populations. In the next section 
we will find it useful to consider a variation to this model, 
which we call the dynamic model, in order to explore the 
sensitivity of our results to the input physics. This change 
is applied to Eq. |2] and assumes that ga s disks are more cen¬ 
trally concentrated at higher redshift ( Mo et alJll998ll and 
are thus more efficient at feeding the black hole during a 
merger. We incorporate this idea in the simplest possible 
way through a transformation of the feeding efficiency coef¬ 
ficient: 0.03 —> 0.01 (1 + z ) (note that the change in coeffi¬ 
cient renormalizes our result to remain on the observed local 
meH - wibuige relation). We point out that now, by construc¬ 
tion, we have introduced an evolution to the 771bh - mbuige 
relation, and this evolution is will be dependent on an evolv¬ 
ing gas fraction. Using both static and dynamic implemen¬ 
tations of black hole growth our interest is to measure the 
strength of the change in the ttibh - mbui ge relation under the 
above physical assumptions to gain a sense for the range of 
possibilities that may seen observationally. 

2.2 A disclaimer 

To keep our model as simple as possible we do not con¬ 
sider other processes through which bulges or black holes can 
grow. This includes the tracking of disk instab i lities which 
contribute to the bulge, as used bv ICroton et alJ ( 2006 1 (this 
is our only variation from their model). Importantly, we do 
not claim that other growth modes are not important to 
the ?tibh - mb u i ge relation. Instead, we assume that mergers, 
as described above, are the primary mechanism that deter¬ 
mines the mass history of the bulge and black hole compo¬ 
nents of a galaxy. This allows us to explore the degree to 
which merger triggered disrupted stellar disks are able to 
drive evolution in the ®bh - mbuige relation. We do not rule 
out the possibility that more complicated processes cancel 
out this effect. 


3 RESULTS 

In Fig. 0 we plot the dibh m-buige relation of our semi- 
analytic model galaxy population at four epochs, z = 
0,1, 3, 6. The filled circles in each panel represent the static 
model described in Section urn while the open squares 
show galaxies where evolution in the black hole feeding ef¬ 
ficiency has been assumed, the dynamic model described in 
Section 12.1.21 For reference, the solid line gives the best fit 


through the observations of iHaring fc Rix 1 (12004 ) for a sam¬ 
ple of 30 galaxies in the nearby universe with well measured 
bulge and super-massive black hole masses. The local popu- 
lation o f both models have been normalized to that found by 
IHaring fc Rixi and thus match it reasonably well. Contrast¬ 
ing this to the z = 6 galaxy population we find clear evolution 
in the m bh - mbuige relation, where even for the static case 
a change in amplitude is observed. For both both static and 
dynamic models this indicates that the characteristic mass 
of a black hole residing in a bulge of given mass is larger 
at high redshift than at low, with the difference between 
models coming only from the degree of evolution found. To 
quantify this evolution we perform a simple \ 2 power-law fit 
with unity slope to each result at each epoch, which we illus¬ 
trate in each panel with long and short dashed lines for the 
static and dynamic models, respectively. These fits are lim¬ 
ited to galaxies with mbuige > 10 9 Mq. We emphasize that, 
for both dynamic and static models, the amplitude of the 
m bh - mbuige relation decreases with time, finally settling 
on the observed relation by the present day. This justifies 
the high normalization chosen in Eq. E Interestingly, both 
models show no significant change in slope with redshift. 

To understand the origin of the evolution found in Fig.E 
we separate the mass growth to model black holes and 
bulges into their respective components. This is measured 
at each redshift by first independently summing the total 
mass contributed from each growth channel (i.e. starbursts 
or disrupted disks for bulge growth, merger driven accretion 
for black hole growth) to all galaxies with bulges having 
mbuige > 10 9 Mq. We then consider the ratio of total bulge 
to black hole growth from these channels (which is also the 
mean relative growth rate) to quantify which dominates and 
when. We do this first for the simplest case, the static model, 
now shown in Fig. E I 11 the top panel we plot the ratio of 
growth rates, mBH/mbuige, for bulge growth through either 
starbursts (dashed line, Eq.E or disrupted disks (solid line), 
both relative to the single black hole growth mode of gas ac¬ 
cretion (Eq.E- (Note that we do not concern ourselves with 
growth from black hole-black hole or bulge-bulge merging, 
since, by definition, this does not change the ijibh mbuige 
relation.) We find that the relative growth of bulges and 
black holes from the existing cold gas supply present dur¬ 
ing the merger is approximately constant with time (dashed 
line). This is expected (see Section 12.1.11 and simply re¬ 
flects the fact that, although individually their growth rates 
can be strongly affected by an evolving cold gas fraction, 
when taken as a ratio this evolution cancels. In contrast, the 
m-BH - mbuige growth ratio from disrupted disks is a strong 
function of redshift, as demonstrated by the solid line, with 
a change of almost an order-of-magnitude between z= 9 and 
the present day. This increase is driven by bulge growth that 
arises from both merged satellite disks and major merger 
disruption of central galaxy disks. 

The bottom panel in Fig. E illustrates how this trans¬ 
lates into an evolution in the msH mbui ge relation itself. 
Here we show the previous \ 2 power-law fits from Fig.E as 
a function of redshift (solid line) and 1 <t scatter around the 
mean (dashed lines). The clear decrease in the amplitude 
of the mBH-mbuige relation by a factor of approximately 
3 tracks closely the increase in the growth of bulges from 
disrupted disks. This demonstrates the simple idea we set 
out in Section E that if mergers are the primary mechanism 
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Figure 1 . The ttibh ~ ^-bulge relation for model galaxies at four different epochs, 2 = 0,1,3, 6. Two realizations of the input physics 
are shown (see Section ED, a “stati c” mo del (circles) and an evolving “dynamic” model (squares). The best fit to the local universe 
observational result of lHaring fc 11 1 >] i2004i) is given by the solid line, while the long and short dashed lines show a simple y 2 power-law 
fit with unity slope to the static and dynamic model results, respectively. These fits have been restricted to galaxies with bulge masses 
m bulge > 10 9 Mq. They highlight a clear evolution in the ttibh m lmlge relation in both models. 


which shape both black hole and bulge growth in the galaxy 
population, then a larger fractional contribution to the bulge 
from disrupted disks should result in an inevitable evolution 
of the ?tibh - rribuige relation. This holds true even when no 
explicit evolution the the growth modes of black holes and 
bulges is assumed. 


In Fig.|!J]we redo the analysis of Figl^Jnow using the dy¬ 
namical model described in Section [ 2 . 1.21 This allows us to 
explore the sensitivity of ttzbh — mbulge evolution when some 
evolution in the black hole growth rate has been assumed. 
The top panel of Fig. [J]clearly shows a much stronger effect 
than that in the previous figure, with the relative growth 
of black holes and bulges from both the cold gas reservoir 
(dashed line) and disk disruption (solid line) changing with 
time with an additional factor of approximately 1 + z. This 
results in a significant boost to the previous evolution seen 
in the ms h -nibuige relation, as shown in the bottom panel, 
with approximately an order-of-magnitude difference now 
predicted between high and low redshift. When one restricts 
the comparison to between redshift z = 1 and the present 
day, the difference in amplitude is still a factor of ~2, which 
may statistically be an observable quantity in the near fu¬ 
ture. 


4 DISCUSSION 

The rapid increase of bulge growth at late times in our static 
model is a consequence of two well understood effects. The 
first is the steady rise in the star formation rate density of 
the universe from high redshift to approximately z = 1 — 2. 
If one accepts, as a general rule, the conventional wisdom 
that the bulk of this star formation occurs in stellar disks, 
then the outcome is a strongly increasing growth of disk 
mass across the galaxy population with time. As disks grow 
the second effect then becomes increasingly important. This 
effect stems from the hierarchical nature of a CDM universe, 
where mergers become more frequent as the universe ages, 
assembling structure from the bottom up. As we discussed in 
Section f2.1 1 mergers also transform disks into bulges. Thus, 
at late times, a larger fraction of the total stellar mass in 
the universe becomes locked up in the spheroid component 
of the population relative to earlier epochs. This results in 
the accelerated bulge growth seen in Fig.|3 and which drives 
the evolution in the 771 bh -mb u i ge relation shown in Fig. 0 
Observationally it is difficult to me asure black hole 
and b ulge masses. In the local universe iMagorrian et all 
(1998) estimate m,BH ~ 0.006 mb u i ge from a sa mple 32 galax¬ 
ies, w hile both iMarconi fc Hunt] |2003) and lHaring fc Rixl 
6 2004 ) independently find wibh ~ 0.002 mb u i ge from im¬ 
proved measurements of ~30 galaxies. Although the statis¬ 
tics are still poor and the uncertainty large, locally at 
least all observations appear to be converging to a con- 
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Figure 2. The static model described in Section 12.1.11 (top) 
The evolution of the relative growth rates for black holes and 
bulges, 772J 3 H /h^bulge? f° r all galaxies hosting black holes with 
bulge masses 772b u ige > 10 9 Mq. As discussed in the text, black 
holes grow from merger triggered cold gas accretion. Bulges, on 
the other hand, grow from both merger induced starbursts and 
disrupted disks, so we plot these two growth channels indepen¬ 
dently (dashed and solid lines respectively). Disrupted disks are 
found to contribute a larger faction to the bulge at low redshift rel¬ 
ative to high, (bottom) The accelerated contribution of disrupted 
disks drive an evolution in the ttibh - 772b u ige relation, shown here 
using the x 2 model fits from Fig.Qacross the entire redshift range. 
The dotted bounding lines show the la scatter of galaxies along 
the relation. 

sistent result. At hi gher redshifts t he pi cture is much less 
clear. For example. Ishields et alJ (120031) claim little evo¬ 
lution in the relation can be inferred out to 2 ~ 3 1 , 
and lAdelberger &; Steidell d2005h measure the quasar-galaxy 
cross-correlation function and find consistency with the lo¬ 
cal ?72 bh - rnbuige ra tio from a sample of 79 2.5 quasars. 

On the other hand, iTreu et "all | |2004|) see variations in the 
ttibh-o' relation at 2 : = 0.37, while iMcLure et alJ d2005h mea¬ 
sure some evolution in the rap h - mbuige r elation using the 
3CRR sample of radio galaxies. Similarly, iRix et al i (GOO l) 
use gravitational lensing to find that quasar host galaxies at 
2 ~ 2 are much fainter than their low redshift cou nterparts 
contai ning quasars of similar luminosity, and [Walter et alJ 
1: 20041 ) find a significant deviation from the local tobh 
w-buige relation for a 2 = 6.4 quasar host galaxy. Future 
observations will need to clarify the exact nature of both 
the high redshift black hole and host galaxy populations. 

Recent theoretical work to understand the cosmolog¬ 
ical assembly of stars an d super-massive black holes have 
led to interesting results. IWvithe fc Loeb ( 12003 ) present a 
model for super-massive black hole growth that successfully 
matches many local and high redshift AGN related observa¬ 
tions. Their work results in a 771 bh — relation constant with 
redshift while predicting the uibh - nibuiga relation evolves as 

1 However see their most recent work i iShields el. al lt2005l) 


Figure 3. As for Fig. [5] however now showing the result for the 
dynamic model, where an evolution in the black hole feeding rate 
has been assumed (Section 12. 1 .21 . Evolution in the uibh ~ m bulge 
relation is now much stronger than that seen previously. 


m-BH/nibuige oc ^( 2 ) 1 ^ 2 (l+ 2) 3 ^ 2 ~ (1 +2) 1 ' 15 , where the final 
approximation is valid when 2 < 2, and £(2) depends only 
on the cosmological parameters and is a weak function of 
redshift. An evolution of this kind would be consistent with 
the scaling assu med in o ur d y namic al model (Section l2. 1.211 . 

Similarly, iMerloni et a.lJ (120041 ) constrain phenomeno¬ 
logically the joint evolution of super massive black holes and 
their host spheroids by fitting simultaneously the total stel¬ 
lar mass and star formation rate densities as a function of 
redshift, as well as the hard X-ray selected quasar luminos¬ 
ity function. With the latter they assume that black holes 
grow exclusively through accretion. Assuming a prese nt day 
disc to spheroid ratio of 0.5 llTasca fc Whitel I 2 OO 5 I) . their 
work favors a model in which the tobh -m.b u i ge relation 
ev olves as ~ (1 + z) 1 ^ 2 . This is a weaker effect than found 
by IWvithe fc Loebl . however demonstrates both the range 
of evolution that may be expected, and most importantly, 
that such non-zero evolution can arise naturally from simple 
studies of black hole and bulge growth. 

As discussed in Section 0 in a ACDM universe the ef¬ 
fect described in this paper will be present in any model of 
black hole and bulge assembly driven by mergers. Indeed, 
this has alread y been seen in the semi-analytic model of 
ICattaneo et alJ |§Q05j) who find similar man -nibuige evo¬ 
lution to that found here (compare their Fig. 6 with our 
Fig- ID • Unfortunately they do not discuss the origin of this 
behavior, but instead choose to focus on the disruption of 
galacti c discs i n relation to the scatter and slope of the re¬ 
lation. ICattaneo et alJ grow bulges both as we do and from 
disk instabilities, which interestingly produces a bi-modal 
utbh - nibuige distribution at high redshift. For simplicity 
we have removed bul ge grow th thr ough disk instabilities 
(as originally used in iCroton et alJ 12006 ). although when 
included we also see such bi-modality. This bi-modal pre¬ 
diction of the high redshift mBH-m.b u i ge relation provides 
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a novel test of the mechanisms through which bulge growth 
may occur. 

Theoretical arguments and numerical work have demon¬ 
strated that galaxy mergers are capable of simultaneously 
triggering growth in both bulges and black holes in a way 
so as to jointly reproduce many of their properties currently 
observed in the local universe. If mergers are the primary 
drivers of black hole and bulge growth in the galaxy pop¬ 
ulation, then we have shown one should expect to see an 
evolution in the tobh - mb u i g e relation which arises from 
an increasing contribution of disrupted disks to bulges as 
the universe ages. In this picture, evolution in the growth 
of bulges drives an evolution in the dibh mbjige relation, 
distinct from the growth rate of black holes. At the very 
least, even if the physics governing bulge and black hole 
growth turns out to be much more complex and cannot be 
expressed in a simplified manner (as is currently assumed 
by most models of galaxy formation), in a ACDM universe 
this effect should still be present and must be included in 
any interpretation of the msH - mbuige relation measured at 
different redshifts. We await future high redshift observa¬ 
tions, e.g. the Gala xy Evo l ution fro m Morphological Stud¬ 
ies (GEMS) project fcix et al.l2fl04l . to clarify the situation 
further. 
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